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Abstract 

We find constraints on heavy, stable, strongly interacting massive parti- 
cles (X) from searches for anomalous nuclei containing them, formed during 
primordial nucleosynthesis. Using existing data, we obtain a limit on the 
abundance ratio Cx = nx/ns in the range of 3 x 10 -8 to 3 x 10~ 13 for masses 
up to 10 TeV if the X — N interaction is sufficiently strong to bind in low Z 
nuclei. We also find a rough lower limit on the X — N interaction that implies 
binding in nuclei with A > 200 over much of the Mx range of interest, and 
address the relative abundance of such anomalous nuclei on Earth. 



Recently a number of authors have entertained the possibility that strongly 
interacting massive particles (SIMPS), with masses a few GeV and above, may 
play a role in particle physics and cosmology |2|, |3|, |], ||, ||, |7|]. For example, these 
particles have been proposed as the source of ultra high energy cosmic rays (UHECR) 
by Chung et al. and others^ H]. There are interesting models of supersymmetry 
breaking where such particles appear as part of the messenger sector or as the gluino 
LSPl. 

Due to their strong interactions and high mass, these particles have distinct 
consequences in cosmology and astrophysics. In order to study their implications, it 
is essential to know their annihilation cross section , a xx which in turn determines 
their relic density It has recently been noted || that, if we assume a xx ~ M x 2 , 
they can appear in sufficient abundance to saturate either the cosmic density or the 
galactic halo. Simple estimates of their relic density give 

Vx = -± ~ 5 x 10 3 -/ 1 

n-y Mpi 

For a 1 TeV X particle, this gives the relative abundance C x e ^ ~ 5 x or 
n x ~ 5 x 10~ 10 / cm 3 . For it to saturate the cosmic density, its mass must satisfy 
M x > 10 3 - 5 GeV. 

If in the process of cosmological evolution, the X particles formed the grav- 
itational potential well and thus became the halo dark matter, their halo number 
density would be considerably enhanced and become n x ~ 3 x 10 _4 /cm 3 or so and 
known estimates of their density could be used to constrain their masses as functions 
of cross sections |^]. 

In [||], using the above form for the X annihilation cross section as well as 

a simple QCD inspired ansatz for the X — A cross it was shown that if the X's 

constituted the halo dark matter and their masses were above 100 GeV, they would 

be in conflict with the results from underground detectors looking for dark matter. 
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The basic reason for this conclusion is that heavy nonrelativistic particles lose only a 
small fraction of their energy in elastic scattering and will therefore have to undergo 
a large number of collisions before they can slow down and get captured. Thus they 
should easily penetrate to reach the underground detectors but their signals have not 
been seen. There are also constraints j5] from the cosmological grand unified photon 
spectrum|| arising from particle-anti-particle annihilation if one assumes either both 
appear in equal abundance or the X particles are their own anti particles. Again 
this would rule out mx > 100 GeV if one assumes the abundance to be halo density. 
Similar conclusions were also derived using the limits on anomalous heavy isotope 
abundances. 

There have been several experimental searches |9|, |Hj to see whether these 
particles exist as dark matter and stringent limits have been placed on their masses 
and their interaction cross sections with matter. Other information on the possible 
existence of SIMPs used in both References |l| and comes from terrestrial searches 
for anomalous heavy nuclei such as those in Ref . . Specifically, the experiment of 
Hemmick et al JTTJ set stringent upper limits on the abundance of nuclei containing 
such heavy stable particles. These limits were used in Ref.'sjl], [| to exclude SIMPs 
above a mass of 1-10 TeV. In order to study this question, one needs to know the 
strength of X — N binding potential Vxn-, which is apriori unknown. One may either 
leave it as a free parameter as we do in the first part or use a "factorization" ansatz 
i.e. a 2 XN = f3a X x a NN which, using the fact that cross sections scale as the square 
of the potential can give a rough idea about the magnitude of Vxn (as we do in the 
second part). 

In this brief note, we present two new results: (A) the first uses an earlier 
investigation of primordial nucleosynthesis of (low Z) anomalous nuclei containing 



X [E3] for the case of sufficiently strong X — N interaction to place limits on the 



X abundance for the case when the annihilation cross section is assumed as in 
Ref. and Vxn is left arbitrary and (B) the second uses the above mentioned 
"factorization" ansatz for the annihilation cross section to determine the likelihood 
of relic Xs binding in high Z nuclei for weaker X — N interaction (Binding is more 
likely for high Z nuclei since their larger radii imply lower kinetic energies.). 
A. Binding in low Z nuclei 

In Reference |L2|, Dicus and Teplitz calculated the abundance of anomalous 



Z > 1 nuclei that would be produced by primordial nucleosynthesis if there existed 
a new neutral, stable, massive baryon X° ( the name used in |12|] was H). This 



was motivated by earlier work of Dover, Gaisser and Steigman|jl3[ on cosmological 
consequences of such particles. These papers assumed that the lightest color sin- 
glet bound state involving the massive colored particle and/or quarks and gluons 
interacts with nucleons in a similar way to that of the A hyperon 

The result of Reference [0 was an estimate of the abundances of anomalous 
nuclei relative to the total present abundance of all isotopes of He, Li, Be, and B for 
the case of low Z binding. It showed a significant enhancement in the abundance of 
9 Be*. Although s Be is not stable, the A hyperfragment 9 Be\ ( 8 Be to which a A has 



been added) is stable. Reference |L2| showed that replacing A by another strongly 



interacting particle, such as X°, present with a cosmic abundance Cx = ^f" = ^f" 
at cosmic nucleosynthesis, would lead to a value for the ratio of the abundance 
9 Be* (anomalous 9 Be) to the sum of all non-anomalous Be isotopes observed today 
greater than the ratio of the abundances of 5 He* to 4 He by a factor of about 10 4 . 

We make the natural assumption that any primordially synthesized anomalous 
heavy isotope should be present in random samples on earth in the same abundance 
as at the time of nucleosynthesis. However to be conservative, we will assume that 
90% of any primordial 9 Be* will have been destroyed in stars. This is a conservative 



estimate because less than 90% of the Universe's Hydrogen has been cycled through 
stars; indeed interstellar deuterium (deuterium is also astrated) is of the order of 
one third the value of primordial calculations. 

We consider (i) the question of what X — Nucleus (also denoted by X — A) 
cross section we expect on the basis of our knowledge of the A hyperfragment and (ii) 
whether for such cross sections, X°'s will be captured during cosmic nucleosynthesis. 



We then apply the constraint of Ref . |L2| to limit either the X — A cross section or 



the X abundance in the light of the data of Ref.s |14] , |15| , \L1 



Following the line of reasoning given for A hyperfragments by Z. Povh||T6l, we 
use for the binding energy of X-particles in nuclei: 



7T~ 



where Vx is the depth (A-independent) of the potential well for the X — A system, 

R = Ro(A — l) 1 / 3 , and u is the reduced mass for the X — A system. We note 

here that X — N binding has been studied in a very helpful paper by R. PlagaJT7|l 

for the case of Mx < 1 GeV. In the A case, R® is about 1.5 Fermi and V\ ~ 27 

MeV. Bi^He^) is about 3.1 MeV. Note that because of the absence of the one pion 

exchange force from the A — N potential, there is no A analog to the deuteron { 2 H\) 

although 3 if a is bound. Based on Eq. (2), we would expect 5 He* which is crucial 

to the formation of Be* and other light anomalous nuclei, to exist at the time of 

nucleosynthesis for some other strongly interacting particle X provided the X — N 

interaction is strong enough to produce an effective square well of depth of 10.5 

MeV. Here we have required that the binding energy Bx be greater than 2.2 MeV 

because primordial nucleosynthesis proceeds rapidly after the temperature becomes 

low enough that the high energy tail of the 7 distribution can no longer dissociate 

deuterium. The low energy (triplet) N-N cross section is about 4 barns whereas the 

low energy A — N cross section is about 0.5 barns. Since the cross section goes as the 
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square of the scattering amplitude which is proportional to the effective potential 
in some approximation, this implies if a(XN) is bigger than 0.1 barns, we expect 
to have a bound state. Although we cannot be certain about the modifications 
to Eq.(2) due to large Mx, it works well for hyper fragments and our intuition is 
that as the kinetic energy decreases relative to the potential due to larger Mx the 
approximation should improve. 

We turn next to the question of whether the X particle will be captured by 
A He nuclei to form 5 He* nuclei. The probability of capture is roughly 

P = n 70 r/ B 0.1(T/T ) 3 ^10- 24 a 6 (3) 

where n 7 o is today's CBR density; tjb is the baryon to photon ratio; the factor .1 
takes into account account the Helium fraction relative to the nucleons; T and t are 
the temperature and time at the time of nucleosynthesis; T is the CBR temperature 
today; and the relative velocity v ~ vn e ~ (6TmJ { 1 e ) 1 ^ 2 ~ Ac (c is the velocity 
of light). Finally cr b is the cross section in barns for X to be captured by A He. 
We compute cr& following Blatt and Weiskopf[|ll| on the electric dipole moment 
contribution to n + p — > d + 7 as discussed in Plaga [|TT|| . For X-capture by a nucleus 
of atomic weight A and atomic number Z. we have 



where 7 = y/(2fj,B), k 2 /2M x ^ .1 MeV and k ~ B. Assuming B ~ 2.2 MeV, 
we estimate that 05 ~ 3 x 10~ 8 barns giving P ~ 10 -3 for large Mx (even for 
Mx — Mx, P — 0.1). Note that this result corrects the assumption of Ref. Jl3| 



copied by Ref. ||12|| that P > 1. It shows that cosmic nucleosynthesis will not "hide" 
stable SIMPs in light nuclei so that signals from XX annihilation discussed in Ref. 
[5] and elswhere must be expected if there is no large X — X asymmetry. 
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Let us now consider the limits on the abundance of the X if the X — N inter- 
action is strong enough for formation of 5 He* and 9 Be*. To derive the abundance 
of 9 Be*, we need to know the abundance of 5 He* relative to normal 4 He at the 
time of nucleosynthesis since the reaction responsible for the synthesis of 9 Be* is 
4 He( 5 He*,^) 9 Be*. One can then write 

H2£|„ = (m*) (^) (2°.) ( 5 ) 

n Be \m He *J \ n B J \n Be J 



From Ref.jll, we get 9 Be* to 5 He* ratio to be lO" 6 ; Eq. (3) tells us that the 
second factor in the above equation is 10 _3 Cj\:, and finally the present Be abundance 



is taken fron Reeves fl~9|j who gives ~ 2 x 10 11 . Combining gives the ratio 



Be* I Be ~ 50CV- Other anomalous isotopes would have abundances below C 



x 



according to Ref.[12|]. In order to reach these conclusions, Ref.[|T^] uses known rates 
for analog production and depletion reactions. 



The data on Be is summarized by Hemmick et al.[n] as follows: Klein et al 
15H have ruled out isotopes of Be up to 93 AMU at concentrations above 10~ 12 



per nucleon using accelerator mass spectrometry whereas Hemmick et al limit the 

abundance to 10~ n or less per nucleon for 100 to 1000 AMU and to 10 -9 or less up 

to 10,000 AMU Q. It should also be noted that Vanderschiff et al. pOfl have placed 

limits on anomalous 5 He* for 42 < Mx < 82. 

Taking into account the above result that Be* /Be ~ 50CV and assuming 90% 

destruction of 9 Be* in stars means that the limits on Cx are: 3 x 10~ 13 to Mx < 93 

GeV, 3xl0~ 12 for 100 < M x /GeV < 1000 and 3xl0~ 9 for 1000 < M x /GeV < 5000 

and 3 x 10~ 8 between 5 < Mx/TeV < 10. These limits will apply to models in 

which a strongly interacting particle X has an attractive interaction with nucleons 

with sufficient strength that the anomalous Be* forms with binding energy greater 

1 ) Profs. Elmore and Hemmick have assured us that the results of Ref.|l|] obtain down to 93 
AMU 



than 2.2 MeV. Roughly such models should include cases in which X — N scattering 
is attractive and <jxn > -1 barns^] • 
B. Binding in high Z nuclei 

We now consider the case in which Vxn is too small for binding with He or Be. 
Even though we do not know the nature of the X — N force from fundamental prin- 
ciples, we can put a lower bound on Vxn using the afore mentioned "factorization" 
hypothesis i.e. a 2 XN = f3a X x a NN which gives 

Vxn ~ Vnn{°'xn / &nn) 1 ^ 2 
> V N N[{a7x l °NN) 1/2 l°NN] 1,2 (3 llA (6) 

For cr™^?, we use the fact that XX annihilation cross-section must be sufficiently 
large to ensure that today Mxrtx < p c - Using standard methodsf|, we find a X -x > 
3 x 1CT 13 barns (independent of M x ) which implies Vxn — 20 KeV if Vnn ~ 50 
MeV. Eq. (2) then gives a condition for the existence of the bound state (in the 
approximation Bx «C Vx)' 

0.425A 2 / 3 /T 1/4 ^ > A- 1 + M x l (7) 
Vnn 

where by Mx we mean M x in GeV's. For A < (M x /GeV), this yields A 5/3 > 

2.3^ / 3 1 / 4 . This yields A ~ 183/? 3 / 20 for R = 1.3 Fermi and ~ 154/? 3 / 20 for 

Ro = 1.5 Fermi. For the case where A 3> (Mx/GeV), we have instead the condition 

for bound state formation A 2/3 > 2] f xU i which yields, A > 427 for M x = 10 GeV. 

Bound states may exist for smaller A and/or Mx if Vxn exceeds its minimum. It is 

important to note that the lower bound on A is very insensitive to the value (3 which 

characterizes any deviation from the formula o\ N = onn&xx- If we se ^ A = 238 in 

2 )Note that if we used the "factorization" formula for cross sections and Vxn ~ 10 MeV, we 
would obtain a much higher annihilation rate in the early universe and consequently a lower relic 
density Cx', nevertheless, it would still be inconflict with the bounds of Ref.Jll|, thus maintaining 
our conclusion. 
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Eq. (8), we should be able to detect bound states of X's with minimal interaction 
down to M x = 83 GeV. 

We can again use Eq. (5) to compute the rate for X-capture. Assuming that 
X-particles have the galactic virial velocity (v ~ lCT 3 c) and halo dark matter density 
(Ph = 5 x 1CT 25 gm cm -3 ), we find for Unranium a ~ 1CT 8 barns. The probability 
of U capturing an X in the lifetime of the solar system is then P = j^uv vir r so i ar ~ 
1.6 x 10- n /(M x (TeV)). 

One can extend the above line of reasoning to discuss what happens when 
the X particles have nuclear cross sections necessary for them to play the role of 
UHECR's. Consider for example the case discussed in Ref. 3 where they need to 
have a o"xat of order of a milli-barn. Using our simple scaling laws, we conclude that 
this implies Vxn — ^ 2 Vnn- Inserting in Eq. (2) gives for large M x that bound 
states should exist for A > 30 (for A <C (M x / GeV)) assuming of course that the 
potential is attractive. Setting A = 238 gives bound states for M x > 6.3 GeV. In 
this case, the capture cross section ranges from 10~ 4 b for M x = 6.3 GeV to 1CT 7 b 
for large M x so that abundances are larger than minimal interaction case (assuming 
similar galactic halo dark matter). 

To date searches for anomalous isotopes only seem to have reached Z = 9 
(Fluorine). Based on the above considerations, we urge experimentalists to search 
for anomalous isotopes with highest Z so that different interaction strengths for 
SIMPs may be explored. In particular, anomalous Uranium may provide a higher 
discovery potential if the halo dark matter is in fact dominated by "minimally" 
strongly interacting massive particles. However searches in much lighter elements 
A>30 should stable SIMPs capable of explaining UHECRs. 
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